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Abstract 

Using all the available empirical information, we analyze the spacing distributions of low-lying 
2 + levels of even-even nuclei. To obtain statistically relevant samples, the nuclei are grouped into 
classes defined by the ratio R4/2 of the exitation energies of the first 4 + and 2 + levels. This ratio 
serves as a measure of collectivity in nuclei. With the help of Bayesian inference, we determine the 
chaoticity parameter for each class. This parameter is found to vary strongly with R4/2 an d takes 
particularly small values in nuclei that have one of the dynamical symmetries of the interacting 
Boson model. 

PACS numbers: 21.60.Ev, 24.60.Lz, 02.50.Tt 
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I. MOTIVATION AND PURPOSE 



During the past decades, a vast amount of nuclear spectroscopic data has been accumu- 
lated. Level schemes involving tens and sometimes hundreds of levels with reliably known 
values of spin and parity are now available for hundreds of nuclei (see Ref. The wealth 
of published spectroscopic data allows for an extensive study of the level statistics of nuclei 
at low excitation energies. In this paper we report on the statistical analysis of low-lying 
states with spin and parity 2 + . 

The interest in such a study derives from the success of random-matrix theory (RMT) 
in describing the spectral properties of nuclear levels (actually: resonances) near neutron 
threshold and proton threshold Careful analysis has shown that the spectral fluc- 

tuation properties of these resonances are in very good agreement with the predictions of 
the Gaussian orthogonal ensemble (GOE) of random matrices. This statement applies, in 
particular, to the nearest-neighbor spacing (NNS) distribution which is well approximated 
by Wigner's surmise jj| 

PwO) = T^exp (-^s 2 ) • (1) 

Here, s is the NNS in units of the mean spacing. In view of the conjecture by Bohigas, 
Giannoni and Schmit 0], the agreement between the spectral fluctuation properties of the 
resonances and the GOE predictions was taken as an indication of chaotic motion in medium- 
weight and heavy nuclei near neutron threshold. Interest then turned to the ground-state 
domain. Here, integrable models often successfully describe the spectroscopic data, and one 
would, therefore, expect the spectral fluctuation properties to be close to those predicted for 
regular systems. For such systems, the NNS distribution is generically given by the Poisson 
distribution, 

p P (s) = exp (s) . (2) 

A statistical analysis requires complete (few or no missing levels) and pure (few or no 
unknown spin-parities) level schemes. Some 15 years ago, complete and pure level schemes 
were available for only a limited number of nuclei (see, e.g., Refs. js> Q|)- The work of 
Ref. Jf| then suggested that the NNS distribution of low-lying nuclear levels lies between 
the Wigner and the Poisson distributions. The evidence presented in Ref. 0] has since 
become an established fact through the work in Refs. js, 10, 11, 12, 13, 14, 15]. 
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The wealth of spectroscopic data now available in the Nuclear Data tables |l( has moti- 
vated us to investigate once again the nuclear ground-state domain. We are able to make 
more definitive and precise statements about regularity versus chaos in this domain than 
has bee. po SSlbl e S0 fa, A s Re, Q, we tocUS atteatioa oa 2+ states o ts e,ec t even-evea 
nuclei. These nuclei are grouped into classes. The classes are defined in terms of the ratio 
R4/2, i-e., the ratio of the excitation energies of the first 4 + and the first 2 + level in each 
nucleus. We argue below that the classes define a grouping of nuclei that have common 
collective behavior. The sequences of 2 + states are unfolded and analyzed with the help of 
Bayesian inference. The chaoticity parameter / defined below is determined for each class. 



II. DATA SET 



The data on low-lying 2 + levels of even-even nuclei are taken from the compilation by 
Tilley et al. Q for mass numbers 16 < A < 20, from that of Endt |l2| for 20 < A < 44, 
and from the Nuclear Data Sheets jlj for heavier nuclei. We considered nuclei for which 
the spin-parity J 71 assignments of at least five consecutive 2 + -levels are unambiguous. In 
cases, where the spin-parity assignments were uncertain and where the most probable value 
appeared in brackets, we accepted this value. We terminated the sequence when we arrived 
at a level with unassigned J n , or when an ambiguous assignment involved a 2 + spin-parity 
among several possibilities, as e.g. J n = (2 + , 4 + ). We made an exception when only one 
such level occurred and was followed by several unambiguously assigned levels containing at 
least two 2 + levels, provided that the ambiguous 2 + level is found in a similar position in 
the spectrum of a neighboring nucleus. However, this situation occurred for less than 5% 
of the levels considered. In this way, we obtained 1306 levels of spin-parity 2 + belonging to 
169 nuclei. The composition of this ensemble is as follows: 5 levels from each of 47 nuclei, 6 
levels from each of 32 nuclei, 7 levels from each of 22 nuclei, 8 levels from each of 22 nuclei, 9 
levels from each of 16 nuclei, 10 levels from each of 14 nuclei, 11 levels from each of 5 nuclei, 
12 levels from each of 2 nuclei, and sequences of 13, 14, 15, 17, 20, 21, 24, 30, and 32 levels, 
each belonging to a single nucleus. 
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III. CLASSIFICATION OF NUCLEI 



A class is defined by choosing an interval within which the ratio 

R 4/2 = E(4+)/E(2+) (3) 

of excitation energies of the first 4 + and the first 2 + excited states, must lie. The width of the 
intervals was taken to be 0.1 when the total number of spacings falling into the corresponding 
class was about 100 or more. Otherwise, the width of the interval was increased. The use 
of the parameter as an indicator of collective dynamics is justified both empirically and 
by theoretical arguments. We recall the arguments in turn. 

(i) Casten et al. Q plotted E(4f) versus E(2f) for all nuclei with 38 < Z < 82 and with 
2.05 < R4/2 < 3.15. The authors found that the data fall on a straight line. This suggests 
that nuclei in this wide range of Z-values behave like anharmonic vibrators with nearly 
constant anharmonicity. As the ratio R4/2 approaches the rotor limit R4/2 = 3.33, the slope 
of the curve showing E(4f) versus E(2f) decreases within a narrow range of £ , (25 f ")-values, 
asymptotically merging the rotor line of slope 3.33. In a subsequent paper [3] it was found 
that a linear relation between E(4f) and E(2f) holds for pre-collective nuclei with R4/2 < 2. 
Thus, from an empirical perspective, the dynamical structure of medium-weight and heavy 
nuclei can be quantified in terms of R\i2- 

(ii) Theoretical calculations based on the IBM-1 model [2(| support the conclusion that 
R4/2 is an appropriate measure for collectivity in nuclei. The model has three dynamical 
symmetries, obtained by constructing the chains of subgroups of the U(6) group that end 
with the angular momentum group SO (3). The symmetries are labeled by the first sub- 
group appearing in the chain which are U(5), SU(3), and 0(6) corresponding, respectively, 
to vibrational, rotational and 7-unstable nuclei. Extensive numerical calculations for the 



classical as well as the quantum-mechanical IBM Hamiltonian by Alhassid et al. |2l( indeed 
showed a considerable reduction of the standard measures of chaoticity when the parameters 
of the IBM model approach one of the three cases just mentioned. The IBM calculation of 
energy levels yields values of R4/2 = 2.00, 3.33, and 2.50 for the dynamical symmetries U(5), 
SU(3), and 0(6), respectively. Thus, we may expect increased regularity of nuclei having 
one of these values of #4/2- 

One might expect that the chaoticity parameter also assumes small values for nuclei 
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near magic numbers, where R4/2 ~ 1- For mass numbers in this domain, our data set is 
unfortunately too small to allow us to draw definitive conclusions, see Fig. ^ 
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FIG. 1: Mean value / (solid lines) and a (error bars) (see Eqs. (|15jl ) of the chaoticity parameter 
for nuclei in several classes defined in terms of -R4/2, obtained by Bayesian inference. 



IV. UNFOLDING 

Every sequence has to be "unfolded", see Ref. j^j], to obtain a new sequence with unit 
mean level spacing. This is done by fitting a theoretical expression to the number N(E) 
of levels below excitation energy E. The expression used here is the constant-temperature 
formula 

iV(£)=iV + exp(^^) . (4) 
We deal with many short sequences of levels. In this case, the unfolding procedure introduces 



a bias towards the GOE. This is shown and discussed in Ref. |23| and will have to be taken 
into account when we discuss our results. The three parameters Nq, Eq and T obtained for 
each nucleus vary considerably with mass number. Nevertheless, all three show a clear ten- 
dency to decrease with increasing mass number. For the effective temperature, for example, 
we find, assuming a power- law dependence, the result T = (15 ± 4)A-(°- 62±0 - 05 ) MeV. This 



result is consistent with an analysis of the level density of nuclei in the same range of excita- 

n 

tion energy carried out by von Egidy et al. [7]. These authors find T = (19 ± 2)y4~( a68±a02 ) 
MeV. 



V. METHOD OF ANALYSIS 

A detailed account of our method has been given in Ref. |23|. Here we confine ourselves 
to the central aspects. We are guided by the idea that the intermediate behavior of the NNS 
distribution of low-lying nuclear levels does not necessarily imply that nuclei in the vicinity 
of the ground state have mixed regular-chaotic dynamics. The key ingredient of our analysis 
is the assumption that the deviation of the NNS distribution of low-lying nuclear levels from 
the GOE statistics is caused by the neglect of possibly existing conserved quantum numbers 
other than energy, spin, and parity. A given sequence S of levels can then be represented 
as a superposition of m independent sequences Sj each having fractional level density fj, 
with j = l,...,m, and with < fj < 1 and Y^JLifj — 1- We assume that the NNS 
distribution pAs) of Sj obeys GOE statistics. The exact NNS distribution p(s) has been 



given in Ref. |2J|. It depends on the (m — 1) parameters fj,j — l,...,m — l. In 25|, this 
expression has been simplified by observing that p(s) is mainly determined by short-range 
level correlations. This reduces the number of parameters to unity and the proposed NNS 
distribution of the spectrum is 

2 

7T e 7T S 

P(s, /) = [!-/ + W)y] exp [-(1 - f)s - Q(/)— ] ■ (5) 

Here, / = Y%=i ff is the mean fractional level density for the superimposed sequences; it is 
the single parameter characterizing the distribution. We determine the function Q(f) from 
the requirement that the expectation value of s is unity, / ds sp(s, f) = 1. This relates Q to 
the error function. We have numerically approximated it and obtain for / in the interval of 
0.1 < / < 0.9 the parabolic relation 

W) = / (0-7 + 0.3/) . (6) 

For a superposition of a large number m of sequences, / is of order 1/m. In the limit of 
m — > oo, p(s,f) — ► p(s,0) = pp(s) as given by Eq. (J2J). This expresses the well-known 
fact that the superposition of very many GOE sequences produces a Poisson distribution. 
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On the other hand, for / — > 1, p(s, f) approaches the Wigner distribution (0) expected for 
a single GOE. We therefore refer to / as to the chaoticity parameter. Our parameteriza- 
tion (|SJ) is not restricted to statistically independent sequences Sj. A system with partially 
broken symmetries can also be approximately represented by a superposition of independent 
sequences In this case, the distribution © which differs from zero at s = 0, is not 

accurate for a domain of very small spacings. The magnitude of this domain depends on 
the ratio of the strength of the symmetry-breaking interaction and the mean level spacing. 

We determine the parameter / by the method of Bayesian inference j^J. Given a sequence 
of spacings s = (s±, S2, sjv), the joint probability distribution p(s\f) of these spacings, 
conditioned by the parameter /, is given by 

N 

p(s\f) = i[p(s l J). (7) 

i=i 

Eq. (J7J) holds if the experimental Si are taken to be statistically independent. This assump- 
tion is justified as long as we confine ourselves to the investigation of the NNS distribution. 
Bayes' theorem then provides the posterior distribution 

of the parameter / given the events s. Here, //(/) is the prior distribution and M(s) = 
1-1 "? (s |/") /if f)df is the normalization. The prior distribution is found from Jeffreys' rule 
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30] 



fJi(f) oc 



1/2 



(9) 



Jp(s\f) [d \np(s\f)/df} 2 ds 

We substitute Eq. (8) into formula Q, evaluate the integral numerically and approximate 
the result by the polynomial 

//(/) = 1.975 - 10.07/ + 48.96/ 2 - 135.6/ 3 + 205.6/ 4 - 158.6/ 5 + 48.63/ 6 . flO) 

Even for only moderately large N, it is useful to write p(s\f) in the form 

p(s|/)=e-W), (11) 

where 

M = (1 - /)<*> + ^ (0-7 + 0.3f)(s 2 } - (ln[l - / + |/(0.7 + 0.3/)a]> . (12) 
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Here the notation (x) = (l/N)J2i=i x i has been used. By calculating the mean values 
(• ■ ■) in Eq. (fT2"j) for various spectra, one finds that the function <j)(f) has a deep minimum, 
say at / = fo- One can therefore represent the numerical results in analytical form by 
parametrizing (p as 

<f>(f) = A + B(f- f ) 2 + C(f - f ) 3 . (13) 

We then obtain 

P(f\s) = cuif) exp(-N[B(f - f ) 2 + C(f - / ) 3 ]) , (14) 

where c = e~ NA /M(s) is a normalization constant. The error interval / ± o 1 ! 2 of the 
chaoticity parameter is defined by the mean value / and the variance a 2 , with 

/ = f 1 fP(f\s) df and a 2 = f\f - J) 2 P(f\s) df . (15) 
Jo Jo 



VI. CHAOTICITY PARAMETER 

The results obtained for / and a are given in Fig. ^ Figure El shows a comparison of the 
spacing distributions conditioned by / and the histograms for each class of nuclei. In view 
of the small number of spacings within each class, the agreement seems satisfactory. 

We recall that the analysis of many short sequences of levels tends to overestimate /. 
Therefore, we focus attention not on the absolute values of / but on the way / changes with 
Ri/2- The graph of / against i? 4 / 2 in Figure 1 has deep minima at R^/ 2 = 2.0, 2.5, and 3.3. 
These values of R4/2 are associated with the dynamical symmetries of the IBM mentioned 
above. Another minimum of statistical significance occurs for 2.25 < R4/2 < 2.35. This 
minimum may indicate that nuclei which lie between the limiting cases of the U(5) and 
0(6) dynamical symmetries, are relatively regular. One may associate this region with the 



critical point of the U(5)-0(6) shape transition in nuclei. Iachello 27( has recently shown 
that this transition is approximately governed by the "critical" E(5i) dynamical symmetry. 
Nuclei with E(5i) dynamical symmetry have R4/2 = 2.2. Experimental examples of this 
critical symmetry have been found by Casten and Zamfir 
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FIG. 2: Comparison of the spacing distributions calculated from Eq. © using the values of / 
given in Fig. ^with the histograms for the empirical NNS distributions for nuclei in several classes 
defined in terms of R^/2- 

VII. SUMMARY 

With the help of a systematic analysis of the NNS distributions for 2 + levels of even-even 
nuclei, we have determined the chaoticity parameter / for nuclei at low excitation energy. 
While in a single nucleus the number of states with reliable spin-parity assignments is not 
sufficient for a meaningful statistical analysis, a combination of sequences of levels taken 
from similar nuclei provides a sufficiently large ensemble. As the measure of similarity we 
have taken the ratio R4/2 of the excitation energies of the lowest 4 + and 2 + levels in each 
nucleus. As seen in Figure 1, the chaoticity parameter / is indeed dependent on R4/2. It 
has deep minima at R4/2 = 2.0, 2.5, and 3.3. These minima correspond, respectively, to the 
£7(5), SO(6), and SU(3) dynamical symmetries of the IBM. A further minimum may relate 
to the critical E(5i) symmetry. 
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